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Abstract. We present evidence that transition region red-shifts are naturally produced in 
episodically heated models where the average volumetric heating scale height lies between 
that of the chromospheric pressure scale height of 200 km and the coronal scale height of 
50 Mm. In order to do so we present results from 3d MHD models spanning the upper con- 
vection zone up to the corona, 15 Mm above the photosphere. Transition region and coronal 
heating in these models is due both the stressing of the magnetic field by photospheric and 
convection 'zone dynamics, but also in some models by the injection of emerging magnetic 
flux. 
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1. Introduction 

Observations of the upper chromosphere and 
of the transition region point to a number of 
mysteries that have yet to be adequately ex- 
plained. 

Chief among these are the rise in the dif- 
ferential emission me asure (DEM) a t temper- 
atures below 10 s K (lMariskal[l992l) and the 
ubiquitous av erage red-shift found in transition 
regio n lines dChae et alJ I1998L iPeter & Judgel 
Il999h In short, the rise in the DEM at lower 
transition region temperatures indicates that 
much more plasma is emitting in this tem- 
perature range than is accounted for by mod- 
Send offprint requests to: V.H. Hansteen 



els where transition region structure is dom- 
inated by thermal conduction from overly- 
ing coronal structures. The average red-shift 
has been explained by a varie ty of models , 
e.g. lAthav & Holzerl d 1982b and lAthavl d 1984b 
conjectured that red-shifts were due the re- 
turn of previously heated sp i cular material, 
on the other hand iHansteenl d 1993b argued 
that nanoflare induced waves generated in the 
corona could produce a net red-shift in transi- 
tion region lines. 

More recentlv lPeter et all d2006) found that 
3d numerical models that span the photosphere 
to corona produce red-shifts in transition re- 
gion lines. In these models coronal heating is 
due the Joule dissipation of currents produced 
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as the magnetic field is stressed and braided by 
photospheric motions, but the u ltimate source 
of the red-shifts is not clear. Spada roet alj 
(2006) used Id models of coronal loops to 
show that red-shifts, and the observed DEM at 
temperatures below 10 5 K could be the result 
of transient heating near the loop footpoints. 

In this work we extend the work of Peter 
et al. and Spadaro et al. to show that transi- 
tion region red-shifts are naturally produced 
in episodically heated models where the av- 
erage volumetric heating scale height lies be- 
tween that of the chromospheric pressure scale 
height of 200 km and the coronal scale height 
of 50 Mm. 

2. 3d simulations of the solar 
atmosphere 

We solve the MHD eq uations using BIFROST 
(Gudi ksen et al.l 1201 Oh . This code employs a 
high-order finite difference scheme and in- 
cludes a high-order artificial viscosity in order 
to maintain numerical stability. These terms 
are also the source of magnetic and viscous 
heating. Non-grey, optically thick radiative 
losses including the effects of scattering are in- 
cluded, as are optically thin radiative losses for 
the upper chromosphere and corona and con- 
duction along the field lines in the corona. 

2.1. Simulation description 

In this paper we will use results from four 
3d simulations of the solar atmosphere. These 
models span the region from 1 .4 Mm below the 
photosphere to 14 Mm above the photosphere, 
and thus include the upper convection zone, 
the photosphere, chromosphere, transition re- 
gion and the (lower) corona. Photospheric and 
lower chromospheric radiative losses are bal- 
anced by an in-flowing heat flux at the bot- 
tom boundary. The chromosphere and regions 
above are maintained by the dissipation of 
acoustic shock waves that are generated in 
the convection zone and photosphere, but per- 
haps more importantly by the dissipation of 
the Poynting flux that is injected into the 
upper atmosphere. This Poynting flux arises 
as a result of the braiding of magnetic field 




Fig. 1. Vertical magnetic field B z in the pho- 
tosphere and magnetic field lines that reach 
the corona (in blue). The image is from time 
t = 2700 s in the A2 simulation. 



lines as well as the injection of new "emerg- 
ing" magnetic flux at the lower boundary. 
These processes have been described earlier 
in the context of numerical simulations of the 
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Three of the models are comprised of 512 x 
512 x 325 points spanning a region of 16.6 x 
16.6 x 15.5 Mm 3 . These models differ in their 
total unsigned magnetic flux, the weak field 
model henceforth named "Al" has (\B Z \) = 
20 G in the photosphere. On the other hand, the 
stronger field model, named "A2" has (| B z |) » 
100 G. In addition, in the stronger field model 
a flux sheet is injected at the bottom bound- 
ary. Finally, the "A3" model is initially iden- 
tical to A2 but with very little flux injected at 
the bottom boundary. The Al and A2 models 
have been run for 1 hour solar time while the 
A3 model for roughly 45 minutes. 

The fourth model, "Bl", has lower resolu- 
tion — 256 x 128 x 160 points — and spans a 
smaller region of 16.6 x 8.3 x 15.5 Mm 3 . The 
unsigned magnetic field is of the same order, 
but slightly weaker (\B Z \) « 75 G, than in as 
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c4 1548.202 A 



Si7 275.350 A 



Line ID 


A [nm] 


log(r max ) 


mass [amu] 


Hell 


25.63 


4.9 


4.0 


C IV 


154.82 


5.0 


12.0 


O VI 


103.19 


5.5 


16.0 


Si VII 


27.54 


5.8 


28.1 


FeXII 


19.51 


6.1 


55.9 


FeXV 


28.42 


6.3 


55.9 



Table 1. List of lines synthesized from the 
simulations described. 



in the strong field model above. This low res- 
olution model also includes an injected emerg- 
ing flux sheet. The Bl model has been run for 
1 hour solar time. 

In figure Q] we show the vertical magnetic 
field in the photosphere along with a subset of 
field lines that reach coronal heights in the A2 
simulation. This magnetic topology is typical 
for the simulations discussed here. 

Also typical for these models is a tem- 
perature structure in the transition region and 
corona where the temperature begins to rise 
some 1.5 Mm above the photosphere and has 
reached 1 .0 MK already at heights of 3 -4 Mm. 
However, the structure of the atmosphere is 
far from horizontally homogeneous andwe find 
very hot plasma (> 2 MK for the A2 run) 
already at heights of 1.5 Mm, while at the 
same time we find cool gas with temperatures 
< 10 kK even at heights up to 6 Mm. 

2.2. Characteristics of simulated TR and 
coronal lines 

Let us now consider some examples of sim- 
ulated spectral lines formed in the transition 
region and corona of these models. We have 
chosen lines that have been previously well 
observed with space borne instruments such 
as SOHO/SUMER or are currently observable 
with Hinode/EIS. These lines cover a wide 
span of temperature in through the outer solar 
atmosphere and are summarized in table Q] 

The synthesized intensities I v in the various 
EUV spectral lines mentioned above are com- 
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Fig. 2. Intensity profiles of the C IV and Si VII 
lines as functions of time from the low resolu- 
tion simulation B 1 . 

puted assuming that the lines are optically thin 
and that the ionization state is in equilibrium 

hv r 

I v = — J 4> v (u, T g )n e n H g(T g )ds (1) 

where the integration is carried out along the 
line of sight, or in this case vertically through 
the box along the z axis. The line profile is 
computed assuming Doppler broadening. 

Examine now how the line profiles vary as 
a function of time at a given typical location 
picked more or less at random in the B 1 simu- 
lation. The transition region line profiles, ex- 
emplified by C IV and Si VII representative 
of the lower and upper transition region re- 
spectively, are shown in figure [2] These lines 
show a fairly similar evolution with significant 
variations in line intensity, shift, and width on 
time-scales down to less than a minute. We see 
strong red-shift events accompanied by bright- 
ening such as at 18 min and at 56 min, and 
a weaker blue-shifted event at 47 min. It is 
clearly evident that the bulk of the C IV line is 
red-shifted at almost all times, while the Si VII 
line has little net shift. 

The coronal line profiles, exemplified by 
Fe XII and Fe XV are shown in figure [3] Also 
here we find large variations in the line pa- 
rameters on several timescales including those 
shorter than 1 min. The coronal lines become 
brighter as the corona grows hotter at times 
later than 20 min in the simulation. The strong 
brightening and blue-shift seen at 47 min is 
found to be associated with a strong heating 
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Fig. 3. Intensity profiles of the Fe XII and 
Fe XV lines as functions of time from the low 
resolution simulation Bl. 



event (the joule heating in the Bl model is 
shown in figure [6] where the event causing the 
line variation is clearly visible). Note that the 
coronal lines at this location are preferentially 
blue-shifted. 

We can also examine the properties of the 
lines by examining the lower order moments; 
the total intensity, average line shift, and av- 
erage line width. In figure [4] we show the line 
moments of the lower transition region C IV 
line at time t = 2700 s in the A2 run. The 
left panel shows the total intensity, the middle 
panel the line shift, and the right panel the line 
width. While the intensities do not obviously 
reflect the topology of the magnetic field, the 
line shift and width are more clearly related to 
the field. Intense red- and blue-shifts are gen- 
erally concentrated to positions where the field 
comes up vertically from the chromosphere be- 
low, i.e in the loop footpoints (see figure [T). 
In the regions between the footpoints we also 
find both red- and blue-shifts aligned along the 
field lines that stretch from the footpoints con- 
centrated in bands near x — 7 Mm and near 
x = 13 Mm. Similarly, we find that the greatest 
line widths are found in the footpoint bands, 
and the C IV lines are in generally much nar- 
rower in the regions between. However, note 
that we also find some large line width regions 
strething in thin strings between the footpoint 
bands and aligned with the field lines shown in 
figure Q] 

The pattern of intensities seen in the coro- 
nal lines is quite different, as seen for exam- 




Fig. 6. Current density squared per particle 
{j]j 2 /p) at time t = 2810 s in simulation Bl 
(red weakest, blue strongest) and vertical mag- 
netic field B z in the photosphere. 



pie in the Fe XII line intensities shown in fig- 
ure[5]which clearly delineate the magnetic field 
lines that stretch more or less horizontally be- 
tween the footpoints in the form of loops. On 
the other hand, the structure of line shifts are 
very similar to the structure found in the tran- 
sition region lines. Blue-shifts are generally 
larger (the maximum blue-shift is -75 km/s 
while the maximum red-shift is only 35 km/s) 
and concentrated to the loop footpoints. Notice 
that footpoint intensities are much lower than 
in the bodies of the loops between the foot- 
points, presumably as footpoint temperatures 
are too low to give emission in coronal lines. 

3. Coronal heating in this model 

Let us now put the "observed" line profiles de- 
scribed into context with the coronal heating 
that characterizes these models. 

The heating per particle (oc j 2 /p) is con- 
centrated to certain regions and structures that, 
in general, follow the magnetic field lines as 
shown in figure [6] where the heating in the 
Bl simulation at time 2810 s. This heating is 
highly irregular both in space and time but is in 
large part concentrated to the transition region 
and lower corona. Note that regions of high 
heating (colored in green — blue) in part are 
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Fig. 4. Total intensity, doppler shift, and line width in the C IV line at time t = 2700 s in the A2 
simulation. The velocity scale is from -40 km/s (blue) to 40 km/s (red). Line widths range from 
narrow black to wide yellow/red with a maximum of 51 km/s. The average line shift in the C IV 
line is 6.6 km/s. 
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Fig. 5. Total intensity, doppler shift, and line width in the Fe XII line at time t = 2700 s in the 
A2 simulation. The velocity scale is from -40 km/s (blue) to 40 km/s (red). Line widths range 
from narrow black to wide yellow/red with a maximum of 76 km/s. The average line shift in the 
Fe XII line is -2.7 km/s. 



aligned in loop like structures, in part in linear 
structures that stretch high into the corona. In 
addition, we find a weaker component (colored 
in red) which, though it also follows the struc- 
ture of the magnetic field, seems much more 
evenly spread out and space filling than the re- 
gions of (sporadic) high heating. In time, the 
linear regions of high heating are short lived, 
of order 100 s or shorter. The linear structures 
sometimes also move horizontally over a thou- 
sand kilometers or so during their lifetime. The 
loop shaped structures are longer lived, though 



the heating varies continuously in time also 
here. 

Thus, high heating events occur in spatially 
localized bursts that lead to increases in the 
temperature but also large velocities. For ex- 
ample, the large linear structure that reaches 
high into the corona, seen in the middle of the 
left panel of figure [6] is the event that is visi- 
ble in figure [2] as a line shift of the transition 
region lines at time t — 47 minutes and that is 
also evident in figure [3] as a line shift and, es- 
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Fig. 7. Average joule heating Tjj 2 /p (solid 
line) and average radiative losses L r = 
-n e nnf(T)/p (dot dashed line) per unit mass 
in the chromosphere and lower transition re- 
gion (upper panel) and corona in simulation 
B 1 . The average temperature is over plotted for 
reference (dashed line). 

pecially for the Fe XV line, by an increase in 
the intensity. 

The concentration of heating towards the 
transition region and lower corona can be 
understood by considering the various scale 
heights involved. We find that joule heat- 
ing caused by the braiding of the magnetic 
field decreases exponentially with height from 
the upper photosphere with a scale height 
that is closely related to the scale height of 
the magnetic energy density (B 2 /2p.o). This 
is the same result as foun d earlier in the 
iGudiksen & Nordlundl (120051) model (see fig- 
ure 7 of that paper), where a heating scale 
height of order 650 km was found in the 
chromosphere, rising to 2.5 Mm in the low 
corona and 5 Mm in the extended corona. 



Similar heating scale heights are found in 
iMartmez-Svkora et all d2008l 120091) . They are 
much greater than the chromospheric pressure 
scale height of 200 km while, on the other 
hand, being much smaller than the transition 
region scale height of 6 Mm or the 60 Mm 
coronal scale height. Thus, in models where 
chromospheric/coronal heating is determined 
by photospheric work on the magnetic field — 
with or without flux emergence — we find that 
the heating per unit mass (or equivalently per 
particle) necessarily is concentrated towards 
the transition region and the low corona: Low 
in the chromosphere the particle density is very 
high and therefore the heating rate per parti- 
cle low. At the same time radiative losses are 
very efficient at these heights, thus the plasma 
easily rids itself of deposited energy. However, 
since the scale height of heating (~ 650 km) 
is greater than that of the particle density (~ 
200 km) in the roughly isothermal chromo- 
sphere, the heating per particle rises exponen- 
tially with height. At some point the heating 
overwhelms the plasma's ability to radiate, and 
the temperature rises to coronal temperatures, 
as illustrated in figures [7] This figure shows 
the joule heating (rjj 2 ) and radiative cooling 
(L r = -n & ertuf(T)) per unit mass as a func- 
tion of height. The temperature continues to 
rise until thermal conduction becomes effective 
enough to balance the deposited energy at tem- 
peratures approaching 1 MK. This temperature 
is high enough that the density scale height 
becomes much larger than the heating scale 
height and the heating per particle therefore de- 
creases exponentially, roughly as the heating 
rate (or the magnetic energy density) itself. 

4. Discussion & Conclusions 

Thus simulations of this sort lead to a coro- 
nal heating model where sporadic heating of 
material occurs mainly in the vicinity of the 
upper chromosphere, transition region, and 
lower corona. When a heating event occurs 
in the upper chromosphere or transition re- 
gion cooler material is rapidly heated to coro- 
nal temperatures. When heating occurs in the 
lower corona the thermal conduction increases 
rapidly with rising coronal temperatures and 
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Fig. 8. Average line shift in the Bl model as 
a function of time for the lines of table Q] (up- 
per panel): He II (dotted), C IV (solid), O VI 
(dash), Si VII (dash dot), and Fe XII (dash dot 
dot dot). The lower panel shows the average 
line shifts and the average velocity (dashed) as 
a function of log temperature in simulation B 1 . 
The error bars in the average line shift are given 
by the standard deviation in the average veloc- 
ities as a function of time. 

the cooler chromospheric material just below 
can be heated by a thermal wave. It is our con- 
tention that these processes happen sufficiently 
rapidly that material is heated more or less "in 
place" before it has time to move. Thus, newly 
heated transition material will be characterized 
by higher than hydrostatic pressures, and we 
expect, on average, to find a plug of high pres- 
sure material that seeks to expand into regions 
of lower pressure on either side; downwards to- 
wards the chromosphere and upwards towards 
the corona. 

The flows away from the high pressure 
plug are reflected in the average velocities 



found in the transition region and lower corona 
as shown in figure [8] Average line shifts found 
in simulation Bl are relatively constant, show- 
ing a large apparent downflow in the transition 
region and a smaller apparent upflow at coro- 
nal temperatures. These line shifts represent 
an actual material flow, on average material is 
flowing away from the temperature of some 
7 x 10 5 K with a maximum downward velocity 
of some 5 km/s at 3 x 10 4 K and a maximum 
upward velocity of 2 km/s at 1 MK as shown 
in the lower panel of figure [8] The line shifts 
reflect this flow trend fairly accurately but due 
to the weighting of the line profile with density 
(oc « e n H the line shifts are exaggerated in the 
transition region. We find that the average line 
shift in the coolest transition region lines, He II 
and C IV are of order 7-8 km/s. The O VI 
line is also shifted, but less, with an average 
that lies near 4 km/s. The upper transition re- 
gion Si VII line is initially blue-shifted, but as 
time passes in the simulation the line shift de- 
creases towards zero and is slightly red-shifted 
at the end of the run. Finally the Fe XII line 
has an average line shift of zero initially but is 
slightly blue-shifted at the end of the simula- 
tion run. 

Note that though a net red-shift is observed 
in transition region lines which corresponds to 
a material flow out of the transition region, the 
coronal mass is actually constant or growing in 
these simulations. Thus, material does not^?ow 
into the corona through the transition region. 
Rather, more material is heated rapidly to high 
temperatures than can be supported, and this 
material is subsequently pushed out of the up- 
per transition region into the corona and chro- 
mosphere. 

Further support for this scenario dominat- 
ing the computed models is shown in fig- 
ure [9] where we plot the correlation between 
the red-shift measured in the C IV transition 
region line and the maximum downward ori- 
ented pressure gradient {dp/dz)/p. The figure 
shows that high red-shift occurs in regions with 
a large downward directed pressure gradient. 
Note that figure[9]also shows that a fairly large 
population of large downward directed pres- 
sure gradients co-exist with line shifts that are 
close to zero. This population is representative 
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Fig. 9. Doppler shift of the C IVline versus the 
maximum downward oriented pressure gradi- 
ent (dp/dz)/p along the atmospheric column 
where the Doppler shift was measured for the 
entire model run of simulation B 1 (see text for 
details). 



of points where the magnetic field is mainly 
horizontal; ie where there is a significant con- 
tribution from the Lorentz force, presumably 
located in the body of low lying cooler loops. 
The latter population does not have much im- 
pact on the average line shifts. 

In summary, we find that outflows from 
the upper transition region can be explained 
as a natural consequence of episodically 
heated models where the heating per parti- 
cle is concentrated towards the lower corona. 
Chromospheric material is heated in place to 
coronal temperatures, or rapidly heated coro- 
nal material just above the transition region 
leads to a thermal conduction front heating 
the chromospheric material just below the 
joule dissipation region to coronal tempera- 
tures; both variations lead to a high pressure 
plug of material at upper transition region 
temperatures that relaxes towards equilibrium 
by expelling material. We have further shown 
that heating with these characteristics naturally 
arises when coronal heating is assumed to be 
caused by the dissipation resulting from the 
braiding of the photospheric magnetic field. 
This mechanism leads to joule heating that has 
an exponential decrease with height closely re- 
lated to the scale height in the magnetic energy 
density. 
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